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ABSTRACT 

We determine star formation rates (SFRs) in a sample of color-selected, star-forming (sBzK) galaxies 
(Kab < 21.8) in the Extended Chandra Deep Field - South. To identify and avoid active galactic 
nuclei, we use X-ray, IRAC color, and IR/radio flux ratio selection methods. Photometric redshift- 
binned, average flux densities are measured with stacking analyses in Spitzer-MIPS IR, BLAST and 
APEX/LABOCA submillimeter, VLA and GMRT radio and Chandra X-ray data. We include averages 
of aperture fluxes in MUSYC UBVRIz' JHK images to determine UV-through-radio spectral energy 
distributions. We determine the total IR luminosities and compare SFR calibrations from FIR, 24 
/Ltm, UV, radio and X-ray wavebands. We find consistency with our best estimator, SFR/fl + [/y, to 
within errors for the preferred radio SFR calibration. Our results imply that 24 fim. only and X-ray 
SFR estimates should be applied to high redshift galaxies with caution. Average IR luminosities are 
consistent with luminous infrared galaxies. We find SFR/j?+;yy for our stacked sBzKs at median 
redshifts 1.4, 1.8, and 2.2 to be 55±6 (random error), 74±8 and 154±17 M Q yr _1 respectively, with 
additional systematic uncertainty of a factor of ~ 2. 

Subject headings: galaxies: high redshift — galaxies: statistics — infrared: general — submillimeter: 
general — radio continuum: general 



1. INTRODUCTION 

The history of star formation traces the origins of visi- 
ble matter in the universe. Understanding star formation 
across cosmic time will yield insights into diverse areas of 
astronomy from the formation and evolution of galaxies 
to the initial conditions of stellar evolution. The redshift 
range 1 < z < 3 is a key epoch in this history, when 
most of the stars in the universe were born. Galaxies 
in this range are identified with color selection methods. 
Estimating their Star Formation Rates (SFRs) is com- 
plicated by contamination from Active Galactic Nuclei 
(AGN), and differences in luminosity to SFR calibration 
in the various wavebands used 13 
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SFRs are estimated from a wide variety of luminos- 
ity calibrations from X-ray through radio wavebands. 
Broadband UV continuum radiation directly probes the 
light of young stars, but is strongly attenuated by dust. 
The thermal IR luminosity, hereafter defined as L/^, 
Liji = L(8 — 1000 /im), mea sures the SFR from the re- 
processed dust emission (e.g. IKennicutt] [l998b'). IR and 
uncorrected UV luminosities represent reprocessed and 
un-reprocessed photons and are used to estimate the to- 
tal SFR. 

UV slope based corrections for dust attenuation (e.g. 
IMeurer et al.lll999f ) have been previou sly applied to high 
redshift, star-forming galaxies (e.g . IReddv fc Steidel 
20041: IReddv et al.ll2005l . [2001 . l20lol : IDaddi et al.ll2007bl : 
Magd is et al.ll2010f ). Low redshift analogs of Lyman 
Break Gala xies (LBGs) have s imilar dust attenuation 
corrections (jOverzier et aL||2011| ), although ULIRGs have 
higher dust correction factors (jHowell et al.ll2010t ) and 
other star-forming galaxies have lo wer dust corrections 
than expected (e.g lBuat et al.ll2010l ). and significant scat- 
ter is observed. 

Radio-wave (1.4 GHz) luminosity in star-forming 
galaxies, primarily synchrot ron emission f rom supernova 
remnants, also traces SFR ([Condon! Il992f ). Radio- wave 
SFR calibr ations rely explicitly upon the IR-radio corre- 
lation (e.g. lB^l2^03l : lY"un et al.ll200lT ) or use an implicit 
conve rsion to Ha luminosity to calibrate SFR ([Condon! 

mn). 

X-ray emission in star-forming galaxies arises from low 
mass X-ray binaries (LMXBs), consisting of long lived, 
low mass (M< 1M ) stars with neutron star compan- 
ions, high mass X-ray binaries (HMXBs), consisting of 
short-lived, massive (M> 8M Q ) stars with a neutron 
star co mpanion, and to a lesse r extent, supernova rem- 
nants (Persic & Rcphacli 2002). These last two sources 

distributions (SEDs). 
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of X-ray emission are linked to short-lived, massive stars, 
providing a rationale for X-ray SFR calibrations. How- 
ever, these multiple sources of X-ray emission, as well 
as X-ray obscuration by gas and dust, complicate X-ray 
luminosity to SFR calibrations. In practice, X-ray SFR 
calibrations (e.g. | Ranalli et al.l [2003: Pe rsic et al.l 120041 : 
ILehmer et all 12010ft explicitly rely upon empirical cor- 
relations with IR luminosity and are therefore indirect 
measures of IR luminosity. X-ray emission from AGN 
are major sourc e s of co nt aminat ion . 

Daddi et al. (|2007bl ) find SFR from dust-corrected 
UV, 24 /im, 1.4 GHz, and X-ray calibrations to be ap- 
proximately consistent f or star-formin g (BzK ; see be- 
low) galaxies at z ~ 2. iReddv et ail (|2006f ) compare 
dust-corrected UV, 24 [iva. and X-ray calibrations for a 
spectroscopic sample of z ~ 2 galaxies, and find dust- 
corrected UV SFR to be consistent with 24 fim SFR for 
Bz Ks, but not other ty pes of star-forming galaxies. 

Wuyts et al. (l20ll compare dust-corrected UV, 24 
[im and Herschel PACS derived SED based SFRs to Ha 
based SFRs and find Ha derived SFRs to require extra 
dust correction to agree with SFR[/y+/j?- I n addition, 
other recent literature has illustrated that 24 /im SFRs 
are overestimated for high luminosity sources at z ~ 2 
(lElbaz et al.UgOll INordon et al.lf2012h . 
'^Ta^meHae ^lir ( 2009ft fi nd ra dio SFR based on the 
calibration of lYun et all ([2001) to be consistent with 
dust-corrected UV SFR for BzKs. 

X-ra y SFRs based on the calibration of iRanalli et al.1 
are found to agree with dust-corrected UV SFRs 
for BX and BM galaxi es in the range 1.5 < z < 3.0 
(IReddv fc SteidellT2004T) and with 24 urn SFR ( Reddv 
et al. 12006ft . ~ 

We seek to expand upon previous studies by includ- 
ing a wider range of luminosity estimates, including the 
most recent results from X-ray observations and analy- 
ses, and by understanding the assumptions and sources 
of uncertainty inherent to each SFR calibration. We seek 
to ascertain the extent to which submillimeter data can 
improve IR-based SFR estimates. We compare differ- 
ent radio luminosity to SFR calibrations, and we use the 
most recent X-ray data and luminosity to SFR calibra- 
tions. By comparing these estimates to the total SFR, 
derived from the sum of IR and uncorrected UV lumi- 
nosity, we seek to test their robustness at high redshift. 

We investigate SFRs binned according to photomet- 
ric redshift, with stacking analyses in radio through X- 
ray wavebands. In particular, we use extensive FIR- 
submillimeter d ata, for which we have an impr oved stack- 
ing algorithm (|Kurczvnski fc Gawiserl 120101 ) . We dis- 
cuss observations and data in Section [2j We present 
our stacking methodology in Section [3] Results are pre- 
sented in Section |4j and we discuss comparisons of SFR 
in Section [5] Our conclusions are summarized in Section 
[HI Throughout this paper, magnitudes are measured in 
the AB s ystem unless stated otherwise. We assume a 
iSalpeterl (l955ft IMF from 0.1 to 1 Mg . Another com- 
monly used IMF, that of iKroupal (|2001l ). would change 
the slope of the low mass end of the IMF, and multiply 
the S FRs presented here by factors of 0.58 (e.g. iHopkinsI 
12007ft . We adopt a cosmology with fi A = 0.7, fl = 0.3, 
and hioo = 0.7. 



2. SAMPLE SELECTION 

The BzK color selection criterion has emerged as a suc- 
cessful color based method for identifying galaxies in the 
range 1.4 < z S, 2-5 in a maximally inclusive manner 
(jDaddi et al.ll2004ft . The quantity BzMj^ is defined as 

BzK = (z - K) AB - (B - z)ab (1) 

Actively star-forming galaxies, sBzKs, are found to sat- 
isfy BzK> -0.2, the upper left region in Figure [1] The 
reddening vector in the BzK plane is parallel to the BzK 
line, making this selection unbiased with respect to dust 
content. In this paper, we ignore the reddest galaxies 
in both z-K and B-z, pBzKs, which tend to be old, pas- 
sively evolving stellar systems, and are located in the up- 
per right region of Figure Q] Comparisons and overlaps 
between BzKs and other color-sel ected galaxy ty pes are 
disc ussed in IReddv et all (|2005l ). iGrazian et all (|2007ft . 
and lGreve et al.1 (|2010ft . SFR estimates ty pically range 
from several tens to hundreds Mm vr" 1 (e.g . iDaddi et all 
27Ml27M[Reddy et alll2005l: JDaddi et al.ll2005l. I2007b |: 
Reddv et al.l|2006l:lDunne et alll2009j:lPannella et all l2009l: 
Greve et allboiol : lYoshikawa et alll2010l ). 

Our sample of sBzK galaxies with Ky eoa < 20 ( K ab < 
21.8) comes from the catalog of iBlanc et al.l (|2008ft . and 
is taken fr om the Multiwavelen gth Survey by Yale-Chile 
(MUSYC; IGawiser et all 120061 ) observations of the Ex- 
tended Chandra Deep Field - South (ECDF-S). Photom- 
etry in the UBVRIz' JHK wavebands is obtai ned from 
the -ftT-selected catalog of Ta ylor et al.l (|2009ft , and is 
augmented with dat a from the MUSYC opt ical catalog 
(Rab < 25.3 depth; iCardamone et al.ll2010ft . which in- 
cludes photometry from 32 bands including 18 medium 
band optical filters and Spitzer Infrared Array Camera 
(IRAC) bands. 

We use X-ray luminosity, IRAC colors and qm criteria 
to discriminate SF galaxies from AGN. There are 110 
BzK sources detected in the combined 250 ks ( Virani 
et al. 120061) + 4 Ms ( |Xue et al.ll2011ft Chandra catalogs. 
Of these sources, 61 are in the Chandra Deep Field - 
South (CDF-S), and the remaining 49 are in the ECDF- 
S. X-ray lumi nosity is used to dis tinguish AGN from SF 
galaxies (e.g. iNandra et alll2002ft : 

L(AGN; 2-10 keV) > 10 42 erg s _1 (2) 

There are 107/110 X-ray detected BzKs that meet this 
criterion (three X-ray detected BzKs in the 4 Ms CDF-S 
have Lx < 10 42 erg s _1 and are considered SF galax- 
ies in this analysis). 25 BzKs were identified as AGN 
on the basis of their position in IRAC color-color space 
([Donley et al.ll2012T ). Finally, we computed qm values, 
defined as the ratio of integrated IR flux, FIR (8-1000 
fim; rest frame; W m~ 2 ) to radi o flux density, F1.4 GHz 
(rest frame; W m~ 2 Hz -1 ) (e.g. Ilvison et al.ll2010bft : 

q IR = log(FIR/3.75 x 10 12 Hz/Fx.4 GHz ) (3) 

14 In this paper, we have corrected colors to account for dif- 
ferences between the Bessel B band filter used on the Very L arge 
Telescope (VLT) in the original sample of Daddi ct al. (2004|) and 
the Johnson B band filter on WFI, Suprime-Cam and other instru- 
ments. Not accounting for this correction can lead to an offset of 
0.5 mag toward lower values in (B - z) color and 0.04 mag in higher 
(z - K) color, producing a significant e xcess of "sBzK" ga laxies that 
are in fact low redshift contaminants (Blanc ct al. 2008). 
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Radio loud AGN are discrimi nated from r adio-qu iet 
AGN / SF galaxies according to llvison etaLl (|2010bl ): 

qm{AGN) < 2.0 (4) 

Computing qj^ for BzKs requires sufficient photome- 
try to estimate the rest-frame FIR SED. Observed frame 
SEDs for BzKs with IRAC+MI PS24+VLA detections 
were fit to iCharv fc Elbadl200ll (hereafter CE01) tem- 
plates. Rest-frame radio luminosities were estimated by 
assuming an S v ~ v a , a = —0-8 r adio SED, character- 
istic of SF galaxies (|Condonl fl992i) . Thirty-five sources 
were identified as AGN using the q/_R criterion; 15 of 
these sources are uniquely identified as AGN with this 
method, and the remainder are also identified as AGN 
using other methods. 

A summary of the results of AGN rejection methods 
is given m Table HJ Sources identified as AGN are ex- 
cluded from the remaining analysis. At present it is not 
possible to conclusively identify all AGN including deeply 
obscured sources; however, these methods make the best 
use of available data. Effects of possible residual AGN 
contamination are discussed in Section [5] 

Photometric or spectroscopic redshifts are assigned to 
sBzK galaxies by matching positions in the K band 
source catalog with the MUSYC optical catalog and us- 
ing redshifts from the optical catalog where redshifts in 
both catalogs are available. Photometric redshifts for 
the o ptical catalog were obt ained with the EAzY soft- 
ware (|Brammer et al.1 1 20081 ): the reader is referred to 
iCardamone et al.l (|2010l ) for details of the redshift deter- 
mination. 

Histograms of photometric and spectroscopic redshifts 
are illustrated in Figure [2] These histograms indicate 
that the majority of sBzKs fall in the traditional redshift 
range associated with BzK galaxies, 1.4 < z < 2.5. There 
are significant numbers of sBzKs at lower (20%)and 
higher (4%) redshifts, although BzK selection is less ef- 
ficient in these ranges. Spectroscopic redshifts are taken 
from the literature, for which there are substantial selec- 
tion biases; therefore we do not expect the distribution 
of spectroscopic redshifts to match precisely the distri- 
bution of photometric redshifts. Thus, the trend of spec- 
troscopic redshifts being distributed at somewhat higher 
redshifts than the photometric redshifts, apparent from 
Figure [2] is not an indication of systematic photometric 
redshift error. 

There are 48 galaxies in the range 1.2 < z < 3.3 that 
have spectroscopic and photometric redshifts available. 
Photometric redshift errors, defined as (z p hot— z spe c) / 0-+ 
Zspec), are well described by a Gaussian with mean = 
-0.02 and er=0.09, see Figure [3] Photometric redshift 
errors for these galaxies are plotted versus spectroscopic 
redshift in Figure [3j and compared to the distribution of 
photometric redshift errors for the larger sample of 1285 
K selected galaxies with both redshifts available, la and 
2cr regions are indicated by shading in the figure. Three 
sBzK galaxies have photometric redshifts that are more 
than 2(T outliers; this outlier fraction is consistent with 
Gaussian statistics given our sample size. 

The final sample of 510 star-forming sBzKs with red- 
shift 1.2 < z < 3.2 are binned according to redshift into 
approximately 1 Gyr intervals in cosmic time: 1.2 < z < 
1.5, 1.5 < z < 2.0 and 2.0 < z < 3.2; details are included 
in Table [2] These sets of galaxies are used in analy- 



sis of individual detections and combined non-detections 
(stacking analyses) to determine their aggregate fluxes in 
each waveband. 

The redshift binned sets are illustrated in a BzK dia- 
gram in Figure [T] (B - z) and (z - K) colors for these 
galaxi es are determined f rom flux values given in the cat- 
alog of lBlanc et alj (|2008f ) . Non-detections in the B band 
lead to undetermined (B - z) colors for 4 and 13 galaxies 
in the redshift bins 1.5 < z < 2.0 and 2.0 < z < 3.2, 
respectively. These galaxies are excluded from Figure Q] 

3. MULTIWAVELENGTH ANALYSES AND STACKING 
ANALYSES 

We measure the redshift-bin-averaged multiwavelength 
SEDs of our sample, and use these data to estimate SFRs 
using published broadband calibrations. Most of these 
galaxies are not individually detected in far IR-radio and 
X-ray wavebands; therefore stacking analysis, using the 
K band positional priors, is essential. The same analysis 
method is used for stacking in each IR-radio waveband: 
an ordinary average of fluxes of individual detections and 
the stacking detection yields an aggregate flux estimate 
for each redshift bin. As discussed in the Appendix, a 
weighted average of individual and stacking detections 
may introduce a bias toward dim sources. The same 
set of galaxies are analyzed in all wavebands with the 
exception of X-ray, in which only sources in the CDF-S 
are studied. We do not expect the additional position 
selection criterion of these sources to introduce any bias. 
Furthermore, as discussed in Section 13.41 galaxies in the 
CDF-S will dominate the stacking signal of sources in 
the full ECDF-S due to the greater X-ray exposure in 
CDF-S. 

3.1. IR-Submillimeter 

MIPS 24, 70 /im. Infrared data were obtained from 
the Spitzer Space Telescope Multi Band Imaging Pho- 
tometer (MIPS) 24 fim and 70 /Ltm images that reach 
5ct d epths of 50 uJv and 3 mJv. respectively f Magnelli 
et al. 120091 ) . The K band sBzK positions are compared 
to the 24 /xm (70 /Ltm) catalog, with positions less than 2" 
(4") separation as the criterion for matching to individual 
detection in each MIPS band, respectively. Individual 
detections are removed from the list for stacking, and 
incorporated into the analysis subsequently. In 24 /jm 
data, in redshift bins from 1.2 < z < 1.5, 1.5 < z < 2.0 
and 2.0 < z < 3.2, there are 64, 90, and 61 individual 
detections and 92, 125, and 78 stacked positions, respec- 
tively. Similarly, in 70 /xm data, in these same redshift 
bins, there are 4, 9, and 1 individual detections and 151, 
205, and 136 stacked positions in each of the redshift 
bins (1, 1, and 2 sources, respectively, could not be ana- 
lyzed due to their being on the edge of the 70 /im image) . 
Stacking is performed on a residual image, after remov- 
ing the matched sou rces from the list, as summarized in 
iHuvnh et alj (|2007l ). The stacking algorithm computes 
an inverse variance weighted average of the flux at each 
stack position. We use an ordinary average of the in- 
dividual detections and the stacking detection to yield 
a single combined estimate for the flux of each set of 
galaxies. We show in the Appendix that, although the 
difference between these two approaches is small, the or- 
dinary average is preferred because the weighted average 
can introduce a bias to the combined flux estimate. 
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BLAST 250, 350, 500 /im. Submillimeter data 
at 250, 350, and 500 /im were obtained from the pub- 
lic archive of the Balloon-borne Large Area Space Tele- 
scope (BLAST) survey of the ECDF-S, which reaches la 
depths of 36, 31, and 20 mJy at 250, 350 and 500 yum, 
respectively, in an 8.7 deg 2 wide field and la depths of 
11,9, and 6 mJy at 250 350 and 500 fim in a 0.8 deg 2 
deep field (|Devlin et al.|[2009f ). 

The redshift binned sBzKs are stacked in the pub- 
lic BLAST 'smooth' data (variance-weighted correlation 
between the signal maps and the effective point-spread 
functions). Each pixel in these data products represents 
the maximum likelihood flux density ( Jy) of an isolate d 
point source centered over the pixel (jTruch et al.|[200l ). 

We use an improved submillimeter stacking and de- 
blending algorithm for stacking in 250, 350, 500 and 
870 /jm data that deals effectively with the problem 
of confusion (jKurczynski fc Gawiserl l20Tol ) . Confusion 
severely limits the effectiveness of s tacking in deep sur- 
veys with limited angular resolution (|Condon| [l974; Hogg 
|2001[ ). particularly at far IR-submillimeter wavelengths, 
and causes a bias in stacking results. Deblending cor- 
rects measured fluxes for confusion from these adja- 
cent sources. This stacking and deblending algorithm 
greatly reduces bias in the flux estimate with nearly min- 
imu m variance. For mor e detail s, see Kurczynski & Ga- 
wiser (120101 ) . All galaxies in the MUSYC catalog with 
Kab < 22 are used in the deblending calculations. 

We find stacking detections (defined as SNR > 3) in the 
250 /im data for the redshift bin 1.5 < z < 2.0 (stacking 
detection SNR = 12), and in the 350 /im data for the 
redshift bins 1.2 < z < 1.5 (SNR = 3) and 1.5 < z < 2.0 
(SNR=10), and in the 500 /im data, for the redshift bin 
1.5 < z < 2.0 (SNR=10). See Table [3] for the stacked 
flux densities and errors in each waveband. In the SED 
fits discussed below, the measured fluxes of formal non- 
detections and their appropriate error bars are included 
in the fits. Combining all of the BzK galaxies, without 
regard to redshift binning, yields stacking detections in 
250, 350, and 500 /an data of 3.9±0.4 mJy, 2.5±0.3 mJy, 
and 1.8±0.2 mJy respectively. 

Fluxes reported from stacking 24 /im selected BzK 
galaxies in the same field, with a different stacking al- 
gorithm are larger by about a factor of tw o than those 
values presented here (jMarsden et al.ll2009f ); in addition 
to the difference in selection of the present sample (which 
is K selected and excludes AGN), this discrepancy may 
also possibly be attribute d to lack of deblendin g in these 
previous reported results ([Chary fc Pope 20 13). 

LESS 870 /mi. Submillimeter data at 870 /im 
were obtained from the Large Apex Bolo meter Cam- 
era E CDF-S Submillimeter Survey (LESS; iWeiss et al.1 
2009), which reaches a la depth of approximately 1.2 
mJy beam -1 . The LESS catalog co ntains 126 individ- 
ually detected submillimeter sources (|Weiss et al.i r2009 N ) 
and these data have bee n used previously for stacking 
analyses of BzK galaxies (jGreve et alJl2010l ). 

The redshift binne d sBzKs are stack ed in the beam- 
smoothed, flux map ()Weiss et al.l 120091 ); galaxies in the 
MUSYC K band catalog are used in the deblending cal- 
culations. Individual detections, as identif ied through 
1.4 G Hz and/or MIPS 24 /jm counterparts (|Biggs et al.l 
120111 ). are excluded from the stacking/deblending analy- 



sis and incorporated into the aggregate (stacking + indi- 
vidual detections) flux estimates as discussed above and 
in the Appendix. There were 2, 6, and 1 individual 870 
/xm detections in the 1.2 < z < 1.5, 1.5 < z < 2.0 
and 2.0 < z < 3.2 bins. These individual detections 
contributed 46%, 29%, 9% to the aggregate (individual 
+ stacking) detection, respectively. Stacked flux esti- 
mates are <282 (3cr), 509±80, and 533 ±100 /iJy in the 
1.2 < z < 1.5, 1.5 < z < 2.0 and 2.0 < z < 3.2 bins. 
These stacked flux estimates are combined with individ- 
ual fluxes into the aggregate values indicated in Table 

EI 

IR Luminosity Estimation. In order to estimate 
Lm, we fit o bserved IR-radiq photo metry to template 
libraries from iCharv fc Elbazl (|200lD . We explore sev- 
eral approaches; in each different region of the 
IR-radio spectrum is chosen for template fitting. Fits 
are performed on each redshift bin-averaged spectrum. 
The template rest frame luminosity is converted to an 
observed frame flux distribution at the median redshift 
of the redshift bin. Then the observed frame model flux 
distribution is convolved with each photometric band- 
pass transmission function to generate predicted model 
photometry. The predicted photometry are combined 
with observed photometric fluxes and errors to generate 
a x 2 statistic for each fit. Each template in the library 
is fit in this way, with the smallest x 2 fit chosen as the 
best fit template. For fits that include observations at 
multiple wavelengths, an overall normalization, A, set to 
its analytical best fit value via d\ 2 jdA = 0, is factored 
into the best fit spectrum. 

We explore several additional approaches to estimat- 
ing Lm from the data using CE01 template fits: for each 
redshift bin, we fit (1) the 24 /im and longer wavelength 
data, (2) the long wavelength IR and radio data exclud- 
ing 24 /im, (3) 24 /im only data. Optical/NIR data are 
excluded from the fits because CE01 templates are con- 
sidered to be incomplete for A < 1/iin. For the 24 /jxn 
only fits, there is no free normalization factor, and the 
CE01 template luminosity is used directly to estimate 
Lm. For these single band fits, the variation of \ 2 with 
template index is used as the basis for determining con- 
fidence intervals; 68% error bars are found for the Lm 
estimate, and these errors are propagated into an SFR 
uncertainty. These fits are illustrated in Table [4] and 
Figured] and are discussed in Section 21 

Uncertainties in Lm are determined from 68% confi- 
dence intervals determined from variations in \ 2 with 
normalization. We considered separately the effects of 
error in the redshift. We approximate the error of the 
median redshift by computing the error of the mean of 
individual redshifts. This error of the mean diminishes 
according 1 / y/~N, where N is the number of objects in 
each redshift bin. We computed Li R for our samples at 
the ±1<7 values of the mean redshift and found the results 
to be the same as the actual Lm to within the normaliza- 
tion error. Consequently, the error in the mean redshift 
does not contribute significantly to the overall Lm un- 
certainty. Because the median redshift is more robust to 
the presence of outliers, we consider the error of the mean 
to be an upper bound to the error of the median. We 
also address the effect of redshift and individual galaxy 
SED errors via simulations, as discussed in Section 13.51 
On this basis, we ignore error of the median redshift in 
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subsequent calculations. 

We address the question of whether Lju determined 
from the average flux SED is indicative of the true aver- 
age of individual galaxy luminosities in two ways: (1) us- 
ing observations of the (bright) subset of 24 /jm detected 
sources, and (2) in simulations, for the entire sample in- 
cluding individually non-detected sources. 

The redshift bin average L/# for individually detected 
24 /j,m sources is computed based on individual fits to 
CE01 templates, using only the 24 fim band photome- 
try. Average Lir is also computed for this 24 /im bright 
subset using the same procedure as used in the larger, 
sample of BzKs: flux values are combined to form an 
unweighted average SED. This average SED is then fit 
to CE01 templates to determine L/^. The two Ljr es- 
timates agree to within 10-20% in our redshift bins. In 
Section f3.5| we generalize these results to our entire sam- 
ple using simulations. 

We use the calibration of iKennicuttl (|1998bf ) to con- 
vert Lir to estimated SFR. This calibration is based on 
the st arburst synthesis models of iLeitherer fc Heckmanl 
(|1995D . and it assumes a continuous burst of age 10-100 
Myr, solar abundances, Salpeter IMF, and bolometric 
luminosity arising from dust reradiation. This calibra- 
tion relates the Lm integrated from 8-1000 fim to SFR 
according to 

SFR IR (M Q yr- 1 ) = 4.5 x 10- 44 L Jfl (erg s" 1 ) (5) 

The uncertainty in this relation arises from uncertain- 
ties in estimation of Lm resulting from extrapolation of 
observed fluxes to total, integrated Lm , confounding 
sources of IR emission that are not associated with star 
formation, and the use of a fixed continuous burst model; 
the combined errors in the SFR are attributed as being 
a factor of ~ 2 — 3. This systematic uncertainty in lu- 
minosity to SFR conversion dominates the overall error 
budget in SFR estimates; the impact of this substantial 
systematic uncertainty on other SFR estimates, many of 
which depend indirectly on the L//J-SFR relationship, is 
discussed in Section [5j SFR estimates for our redshift 
binned sBzKs are discussed in Section 14.11 Uncertain- 
ties in these SFR estimates include only errors in Lj# 
(arising from photometry, and template normalization). 

3.2. UV-Optical-NIR 

MUSYC 5cr imaging depths include U = 26.5, B=26.9, 
V=26.6, R=26.3, 1=24.8, z'=24.0, J =23.1, H=22.4, 
K=22.4 as well as 18 medium band photometry in the 

range [4270,8560] A (IGawiser et all [20061: 1 Taylor et all 
2009; ICardamone et alJ 120101 )" Galaxies are individu- 
ally measured in MUSYC UBVRIz'JHK bandpasses 
via aperture photometry, and additionally their fluxes in 
each redshift bin are combined in an unweighted average 
to yield a single averaged SED for each redshift range. 
These average fluxes for redshift binned sBzKs in each 
UV-radio waveband are indicated in Table [3] IR data in 
the ECDF-S are available in IRAC bands at 3.6, 4.5, 5.8 
and 8.0 ^m (SIMPLE: [Damen et aT1l2011l ). As discussed 
in Section [21 IRAC data were used in photometric red- 
shift determination; however, these data were not used 
in the fits to determine Ljr. 

To determine UV continuum luminosity before dust 
correction, L^ ncorr , we use the available optical-NIR 



photometry. We estimate the rest frame 1600 A flux 
density, fUncorr^ y j a interpolation of the two bracketing 

o 

broadband fluxes. The specific luminosity at 1500 A , 
jjjncorr j n un [^ s f e rg s _1 Hz -1 , at the redshift, z, is 

then found from the flux density, f^ ncorr in units of /xJy 
and the luminosity distance, Dl, according to 



jjjncorr ^ x 1Q~ ^ jUncorr ^^^L 



(1 + *) 



(G) 



To convert lumi nosity to SFR, we use the calibration of 
IKennicuttl (Il998af). wh ich corresponds to the calibration 
of IMadau et all ([19981 ) converted to Salpeter IMF and 
0.1 - 100 M mass limits. 

SFR UV (M Q yr' 1 ) = 1.4 x lO^L^erg s^Hz' 1 ) (7) 

This calibration assumes continuous star formation over 
timescales of 10 8 years or longer and solar metallicity. 
IKennicuttl (j!998a[ ) discusses sources of systematic uncer- 
tainty in the various published L v - SFR calibrations as 
arising from the use of different stellar libraries and as- 
sumptions about the star formation time scales; the pub- 
lished calibrations differ by about a factor of 2. We use 
Equation [7] to compute SFR from uncorrected luminosi- 
ties, L^ ncorr , to estimate the contribution to SFR that 
is unobscured by dust. We also apply this expression to 
dust-corrected luminosities to determine dust-corrected 
UV SFR, SFRg^ rr , as discussed below. 

To obtai n dust-corrected UV SFRs, we use the method 
of IRX-g (|Meurer et aTlll999D . T his method has been 
used in high redshift LBGs , (e.g. IMeurer et all 119991 : 
lAdelberger fc Steidell 120001 : IReddv et al.l 120101 ) . BzK 



galaxies (IDaddi et al.lffi wb) , and galaxies at lower red- 
shift (e.g. iBuat et al.ll2010l: iHowell et al.ll2010l: Takeuchi 
et al. 120101) . 

We fit the rest frame UV f\ spectrum to a power 
law, f\ (A) = AX 13 , using a range of trial (3 values: - 
2.5</3<1.0in steps A/3 = 0.01. Due to the availability 
of 18 medium band photometry, these fits typically had 
12-14 sampled points in the spectrum. The wavelength 

o 

range of the fits, [1268, 2580] A in the rest frame, is 
chosen to be the same as t hat u sed in Cal zetti et al.l 
(|1994[ ), which IMeurer et al.l ([1999D also adopted. These 
values are redshifted into the observed frame, and pho- 
tometry data falling within this range are used for the 
fits. For each trial value of (3, predicted flux values are 
computed at each relevant, observed wavelength by in- 
tegrating over the filter bandpass transmission function, 
T(A), and i ntergalactic m edium transmission function, 
M(A), from IMadau! (|1995l ). The integral is expressed in 
terms of the number of photons detected, hence an extra 
factor of A is included in the integrand, as illustrated be- 
low in Equation [8] The integral is normalized to units of 
(ijy by dividing by the corresponding integral of a refer- 
ence spectrum that is flat in /„ (/£ e ^=l fiJy), which is 
converted to a photon number spectrum. This approach 
leads to the expression for predicted flux density, /£T , 
for each broadband filter, i 

fr\,Jv)-A^" T f^ XdX (8) 

iv% ^ J/r / (A)T l (A)AdA 
We compute a % 2 for each fit, and we optimize the nor- 
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malization parameter, A, by selecting the value for which 
dx 2 /dA = 0. Finally the complete, normalized predicted 
flux in fijy is computed from Equation [8] The above fit 
procedure is repeated for each trial /3 value, and the fit 
with the smallest % 2 is chosen to represent the data. 

The resulting power law index, /3, is then used to com- 
pute the UV extinc tion from the empirical relation of 
iMeurer et ail (|1999f ) 



A 1600 = 4.43 + 1.99/3, 



(9) 



which is found to have 0.55 mag dispersion about their 
fit in Ai6oo and a standard error in the fit zero point of 
0.08 mag: see Equation (11) and Figure 1 from Meurer 
et al. ([19991 ) . The UV extinction is then used to correct 
the measured UV flux according to 



F, 



Corr -^q0.4Ai6oo xpUncorr 



UV 



UV 



(10) 

Finally the corrected UV flux is used to estimate the 
SFR using Equation [7] The above procedure is executed 
for each sBzK galaxy individually as well as for the un- 
weighted average spectrum of all sBzKs within a redshift 
bin, whereby the bin median redshift is used to compute 
the luminosity distance. 

Uncertainties to the corrected UV luminosity arise 
from observed flux uncertainty, error in the dust correc- 
tion factor, and photometric redshift error (for individ- 
ual galaxies, we adopt the value of Sz/(1 + z) = 0.009; 
for the average spectra, this contribution is negligible, as 
discussed in Section 13.11 and via simulations in Section 
13. 5[) . These uncertainties are combined using standard 
error analysis. We do not include the systematic un- 
certainty associated with the luminosity— SFR calibra- 
tion (discussed above, about a factor of 2); rather we 
consider this uncertainty with similar systematics from 
other waveband estimators separately in Section [5j The 
results of these computations are shown in Table [5] and 
discussed in Section [5] 

3.3. Radio Luminosity and SFR Estimation 

Radio. Radio data at 610 MHz were obtained from the 
Giant Metrewave Radio Telescope (GMRT) survey of the 
ECDF-S, which reac hes a typical depth of 40 /xJy beam -1 
(jlvison et al.l I20l"0ah . 1.4 GHz data were obtained from 
the Very Large Array (VLA) survey, which covers the 
ECDF-S to a typical d epth of 8 /xJy beam -1 and includes 
464 cataloged sources (jMiller et al.ll2008l ). 

Flux estimates are found from weighted average im- 
age stacking (excluding individual detections, which are 
included after stacking) as well as median image stack- 
ing (of all sources) of the VLA and GMRT data. Me- 
dian stacking is commonly used to reduce the influence of 
radio-loud AGN. We adopt the weighted average method 
to be consistent with the analysis in other wavebands, 
and we adopt median stacking for comparison. In the ra- 
dio regime, where the spatial resolution is relatively high, 
making images allows us to conserve flux density that 
would otherwise be lost due to smearing by astromet- 
ric uncertainties and finite bandwidth (chromatic aber- 
ration) at the cos t of larger flux density uncertainties 
(jlvison et al. 2007). Radio fluxes, luminosities and cor- 
responding SFRs are illustrated in Table [S] 

Our data include flux measurements, S„, at 1.4 GHz 
and 610 MHz for each of the redshift binned SEDs. The 



radio spectra l index, S^ o c v a , (typical a ~ -0.8 for 
galaxies, e.g. lCondonlll99 2f) is estimated for each of our 
three redshift bins to be -0.74±0.2, -1.20±0.2, 0.06±0.4 
respectively. In the higher redshift bins, these computed 
indices deviate significantly from the a = -0.8 for syn- 
chrotron emission (e.g. anomalously high 610 MHz flux 
estimate in the 1.5 < z < 2.0 bin). Similarly high fluxes 
were also reported fo r galaxies in the range < z < 2 in 
iBourne et al.l (|2011[ ). and interpreted as resulting from 
AGN contamination at high redshift. Therefore in keep- 
ing with other reported literature, we adopt the a = -0.8 
value in computing luminosities and SFRs. 

In estimating radio luminosities, we use the median 
redshift, z, and the corresponding luminosity distance, 
Dl, in Mpc to compute the aggregate rest frame 1.4 GHz 
luminosity, L v> xa GHz in units of W Hz -1 , from the ob- 
served frame 1.4 GHz flux, S v , in units of /iJy according 
to 



LuAA GHz = 9.523 x 10 1 



(1 



*i/,1.4 GHz (ttj 



)l+a 



To est i mate SFR from gh ? i we use the model o f 

iCondonl (119921) as imple mented in lHaarsma et all (2000) 
and lDunne et all (|2009l ). Following the implementation 
in lHaarsma et all (|2000D . SFR in units of M yr" 1 is a 
function of frequency in units of GHz, L„ 1.4 GHz in units 
of W Hz^ 1 , scaled by a factor Q and is given by 



C p nC ondon 



Jv,XA GHz 



5.3 x 10 21 i/ Q + 5.5 x 10 20 ^- 01 



(12) 



We use th e value Q = 5.5 to scale the SFR(M>5 M ) cal- 
culated in lCondonl (fl992l ) to SFR(0.1-100 M ) used here; 
this scaling factor de pends on the assumed (Salpeter) 
IMF used here and bv lHaarsma et al.l (|2000l ). 

For comparison, we also e stima t e SFR from 1.4 GHz 
flux using the calibration of iBelll (120031 ) . This calibra- 
tion is based on the IR-radio correlation; it assumes that 
nonthermal radio emission directly tracks the SFR, and is 
chosen so the radio SFR matches the IR SFR for L > L* 
galaxies. The SFR calibration 



SFR?^ Hz (M @ yr' 1 ) = 5.52 x lO" 2 " 



L 



v.lA GHz 



(13) 



is adopted here. A similar calibration is found in Yun 
et al. (120011). SFRV °" a g" exceeds SFRf ^ " „„ by a factor 
of two; the calibration of ICondonl (|1992D explicitly mod- 
els the thermal and non thermal em ission mechanisms, 
whereas the calibration of Bell (2003) relies upon the IR- 
Radio correlation. Thus we expect agreement between 
SFRf 

AGHz an d IR-based SFR estimates, if the IR-radio 
correlation continues to hold at high redshift, as has in- 
deed been suggeste d in the literature (Sa rgent et al.ll2010l : 
llvison et ai1l2010bl ). 

Uncertainties to the radio luminosities are computed 
by incorporating uncertainties from redshift and flux 
measurement; these uncertainties are propagated into 
the SFR uncertainties. When only flux measurement un- 
certainties are included in the error budget, uncertainties 
in SFRf agree to within 30% of published value s 
punne et all [2009). In the calibration of HH3 (|2003j), 
scatter in the IR-radio correlation contributes a factor 
of 1.8 (dispersion of 0.26 dex for individual galaxies) to 
the uncertainty and dominates the total error budget; 
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this additional systematic uncertainty arising from the 
L/fl-SFR calibration is discussed in Section [5] 

3.4. X-Ray SFR Estimation 

X-ray exposure in the ECDF-S consists of 4 Ms in the 
central w 16' x 16' CDF-S, reaching approximate sen- 
sitivities of lxl0~ 17 and 7xl0~ 17 erg cm" 2 s" 1 in the 
0.5-2.0 and 2.0-8 keV bands, respectivel y, and giving thi s 
field the deepest X-ray coverage to date (|Xue et al.ll2011[ ). 
These data are augmented with four flanking 250 ks ex- 
posures that complete the « 30' x 30' ECDF-S field, and 
reach sensitivity limits of 1.7xl0~ 16 and 3.9 xlO -16 erg 
cm" 2 s" 1 in the 0.5-2.0 and 2.0-8.0 keV ban ds, respec- 
tively ((L^^EBISIlffll Iw^f^ISIlSS)- 

X-ray stacking analysis was performed in the 4 Ms 
CDF-S; due to the ratio of exposure times in the CDF-S 
versus ECDF-S, the deeper CDF-S data will dominate 
any stackin g signal. The X-ray stacking algorithm is 
discussed in lTreister et al.l (|2011[ ); a position dependent 
aperture correction was used to account for the vary- 
ing Chandra PSF with off-axis angle, and to minimize 
this correction, only sources within 10' of the aim point 
were stacked. Sources that have an X-ray detection closer 
than 15" to the stacking position are removed to provide 
a better estimation of the background. This procedure 
leaves 19, 29, andl9 source positions in redshift bins from 
1.2 < z < 1.5, 1.5 < z < 2.0 and 2.0 < z < 3.2 respec- 
tively that are stacked in Chandra soft band and hard 
band data. Stacked fluxes are combined with one in- 
dividually detected source each in the 1.2 < z < 1.5 
and 1.5 < z < 2.0 bins, according to the procedure in 
the Appendix. Including these individual detections in- 
creased the soft band stacking flux estimate by 33% in 
the 1.2 < z < 1.5 bin and had negligible effect on the 
1.5 < z < 2.0 bin. 

Counts to flux, and flux to luminosity conversions are 
done assuming a spectrum with photon index T=1.2 and 
cutoff energy, E c =20 keV. The rest frame 2-10 keV lu- 
minosities are used to es timate SFR. As shown by Persic 
& Rephaeli (|2002l ). the X-rav spectrum of star-forming 
galaxies that do not have an AGN is dominated by 
HMXBs, which are best described by T = 1.2 and a 
cutoff energy of 20 keV. Many star-forming galaxies also 
present a thermal componen t, which is typic ally softer in 
X-rays with kT ~ 0.7 keV (|Fabbiandll989T ). The spec- 
trum of the resulting combination is something softer 
than a pure T = 1.2, but n ot quite r = 2. In studying 
LBGs, iNandra et all (|2002l ) assumed an intrinsic spec- 
trum of F = 2.0, more typical of local Seyfert galaxies 
and soft X-ray selected quasars. To estimate the effect of 
different assumptions of photon spectrum index, we also 
computed counts to flux and luminosity conversions us- 
ing r = 1.9. The differences of conversions from counts 
to flux are ^11%. Similarly, the differences in conversion 
from observed frame soft band to rest frame hard band 
are ~ 16% between these two assumptions of spectral in- 
dex. Thus in the soft band, the uncertainties due to an 
assumed spectral shape are ~ 20%. X-ray fluxes, lumi- 
nosities and SFRs are tabulated in Table [7] and discussed 
in Section l4~4l 

The X-ray - SFR calibration of iRanalli et~al1 (120031 ) is 
widely used, and is based upon the X-ray - IR luminosity 
correlation observed in galaxies with L2-iofe e y ^ 10 41 erg 
s" 1 . The SFR in units of M yr" 1 is related to the 2-10 



keV luminosity, L2_iofc e y in units of erg s , according 
to 

SFR^LeV = 2-0 x 10- 40 L 2 _ lofce y. (14) 

This ca libration implicitly assumes the L/ij-SFR calibra- 
tion of lKennicuttl (|l998al ) and a Salpeter IMF (0.1-100 
Mq) consistent with other calibrations mentioned in this 
paper0 However, the Ranalli et alJ (|2003l ) sample in- 
cludes few star-forming galaxies in the ULIRG regime, 
where the Ly - SFR correlation is observed to drop 
(|Lehmeill20l"0f ). 

The uncertainty to SFR^" a ^ is computed by adding 
in quadrature uncertainties in X-ray luminosity and the 
0.09 dex error of the slope in the X -ray— IR- luminosity 
correlation fsee IRanalli et all 120031 Equation 10). Lu- 
minosity uncertainties are computed by propagating the 
flux estimate errors; redshift errors of the average spec- 
trum can be neglected, as discussed in S ectio n 13.11 and 
shown in simulations discussed in Section 13.51 

Subsequent studies have related instantaneous SFR 
specifically to luminosity from short-lived (e.g. Grimm 
et al. 20031: IColbert" et al .1 120041 : IPersic et al.H2004h while 
slowly evolving, LMXBs are linked t o stellar mass, i.e. 
integrated S FR dColbert et al] l20"olh. The X-ray SFR 
calibration of IPersic et a l. (2004) is based upon the lumi- 
nosity of HMXBs, and it relates SFR in units of Mq yr _1 
to the 2-10 keV HMXB luminosity, ^"-iok^v in units of 



erg s , according to 



C TP jjPersic -i n 

or -n-2-lOfceV — iU 



-39 t HMXB 
lj 2-WkeV- 



(15) 

The fraction, /, of HMXB X-ray luminosity to total 
X-ray luminosity has been estimated as / ~ 0.2 (with 
substantial scatte r due to low statist ics) for nearby star- 
forming galaxies (jPersic et alJ 12004). For high redshift 
[z > 1) galaxies, in the absence of definitive estimates 
from X-ray spectroscopy, the value / = 1 has been used 
on the assumption that LMXBs (or other sources of emis- 
sion) contribute a n egligible fraction to the total X-ray 
lumin osity at z ~ 2 (jPersic et~aT1[2004l : rPersic fc Rephaehl 

HmzI). 

Assuming / = 0.2 for nearby star-forming galaxies 
A% = 0- 2 L ?-wk,v, which brings the cal- 
intq equivalence with the 
20031 ). Our data for high 



calibration of IRanalli et al 
redshift galaxies do not directly constrain the HMXB lu- 
minosity fraction; with the assumption / = 1 



leads to Ll 



ibration of IPersic et alJ (1200 



an d the SFRs estimated from IPersic et al.1 



(2004) exceed those of lRanalli et"aLl 1)20031 ) by a factor of 
5. In comp uting SFR, we adop t / = 1 for the X-ray cal- 
ibration of iPerskTe t al. (2004j) for our sample of sBzKs, 
and we regard the resulting SFRs as upper limits. 

The relative contribution of LMXBs to the total X-ray 
lumin osity is believed to decline above z ~ 1 ( Ghosh & 
White |2001|)7 and to be subdominant in high SFR (e.g. 
> 100 M© yr" 1 ) galaxies. A bilinear relation between 
X-ray luminosity and b oth SFR and stella r mass, M*, 
has been proposed, (e.g. iColbert et "aT1l2004D . The X-ray 
SFR calibration of iLehmer et al.l ((20101 ) isderived from 



15 N ote that the radio SFR calibration cited in R analli et al.1 
(2003) refers to M> 5 M Q mass range. For a Salpeter IMF, the 
resulting X-ray-radio derived SFRs differ by a factor of 5.5 from 
the 0.1-100 Mq range used here. 
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such a relationship 

LtioCeV = aM* + PSFR. (16) 

In analysis of LIRGs/ULIRGs exte nding to L 2 -iofc e v 
~ 10 41 - 5 erg s~^ lLehmeret~aII(p0Tl report a = (9.05 ± 
0.37) x 10 28 erg s" 1 Mq 1 and ft = (1.62 ± 0.22) x 10 39 
erg s _1 (M Q yr" 1 ) -1 . In the absence of the M* term 
in Equation [T6l and with the assumption / = 1, as 
discussed above, this calibration b ecomes consi s tent t o 
within errors of the calibration of IPersic et a l. (2004), 
after accounting for the differences in IMF assumed by 
these authors. 

In order to compute SFRjf ^oleV' we estimate stellar 
masses for our redshift binned sBzK samples. We use the 
empirical correlation between the observed frame K band 
magnitude and the stellar mass for sBzKs at z > 1.4, 
determined from SED fits, that is presented in Daddi 
et al. (120041) . 

log(M*/10 u M Q ) = -0A(K tot - K 11 ) (17) 

where K 11 = 21.4 is the K band magnitude correspond- 
ing on average to a stellar mass of 10 11 M Q . IDaddi et al.l 

(2004) report uncertainties of ^40% with this relation. 
Uncertainties in SFR^^Jy are calculated by prop- 
agating the uncertainties in the luminosity and stellar 
mass, along with the reported uncertainties in the pa- 
rameters a and ft indicated above. This cali bration is 
based implicitly on the L/^-SFR calibration of lBell et al.l 

( 2005) which yields lower SFR s by ~13% compa red to 
the corresponding calibration of lKennicutt] (|1998bl ); how- 
ever, we neglect this small calibration difference, and dis- 
cuss systematic uncertainties in comparison with other 
SFR indicators in Section [5] 

3.5. SIMULATIONS 

We investigate the effects of redshift bin averaging, 
photometric redshift errors and dispersion of individual 
galaxy SEDs on our stacked L/^ estimates with simu- 
lations. Averaging the photometric flux density from 
galaxies at slightly different redshifts within a redshift 
bin introduces redshift smearing. In order to study this 
effect, we simulate a set of identical CE01 template spec- 
tral models. These spectra are shifted to the identical 
redshifts of galaxies in our 1.5 < z < 2.0 redshift bin, and 
then averaged together, analogous to the actual stacking 
procedure. The redshift bin averaged spectrum is nearly 
identical to the template spectrum except for moderate 
smearing of the emission peaks that contributed a negli- 
gible amount to the integral; consequently, the quantity 
of interest, the integrated IR luminosity, is robust against 
redshift smearing over our bin widths. 

Of greater concern is the effect on L/^ due to photo- 
metric redshift errors and dispersion due to individual 
galaxy SEDs. To quantify the contribution of these er- 
rors to the estimated L//j for each redshift bin, we sim- 
ulate sets of galaxies with SEDs chosen at random from 
CE01 templates and distributed in redshift to simulate 
the observed source distribution. The photometric red- 
shift error distribution is determined from comparison 
of spectroscopic and photometric redshifts for the subset 
of sources with both estimates available and is shown in 
Figure |3l This distribution is well fit by a Gaussian with 



mean = -0.02 (i.e. bias) and a=0.09 (i.e. scatter). The 
bias is first removed from the simulated object redshifts, 
and then artificial redshift errors drawn from this biased, 
Gaussian distribution are added in each repeated trial of 
the simulation. The resulting spectra are averaged, and 
this averaged spectrum is integrated to determine L/jj. 
An example of these spectra from the 1.5 < z < 2.0 red- 
shift bin simulation is illustrated in Figure [5j From the 
figure, it is apparent that the averaged spectrum has a 
slightly higher flux than the single object spectrum near 
the emission peak and therefore will overestimate Lir. 
For each repeated trial, the fractional error in the 
estimate is determined by comparing the bin averaged 
Ljb. to the true redshift bin averaged luminosity. The 
frequency distribution of Ijir fractional errors is deter- 
mined directly from 10 4 repeated trials for each redshift 
bin. The fractional error distributions for redshift bins 
1.2 < z < 1.5, 1.5 < z < 2.0 and 2.0 < z < 3.2 are each 
Gaussian with mean = 12%, 14%, 2% (bias) and er = 
7%, 6%, 6% (scatter) respectively. The fractional error 
distribution for the 1.5 < z < 2.0 bin is illustrated in 
Figure [5j Our reported values in Table 0] are bias sub- 
tracted, and have the scatter added in quadrature with 
other sources of errors. 

Finally, in simulations we address the issue of whether 
luminosities computed from the average flux SED taken 
to be at the median redshift of each bin may accurately 
reflect the true average luminosity of our sample of indi- 
vidual galaxies. To test this method, we distribute a set 
of galaxies, with CE01 templates chosen at random, dis- 
tributed in redshift according to the actual source pop- 
ulation, compute the of each galaxy individually, 
and average them to determine true average Then 
for each galaxy, we compute observed frame photomet- 
ric fluxes in each FIR-radio waveband, and subsequently 
compute the sample average observed flux SED. We com- 
pute the L/ij of the average flux SED, assumed to be at 
the median redshift using the methods of Section lBTTl We 
compare this estimate with the true average I^ir, 
and determine the distribution of errors with 10 3 Monte 
Carlo realizations. 

The fractional error distributions for redshift bins 
1.2 < z < 1.5, 1.5 < z < 2.0 and 2.0 < z < 3.2 are 
each Gaussian with mean = 0.02, 0.09, and 0.02 respec- 
tively and a — 0.03 in each case. Thus the average flux 
spectrum approximation introduces only a small redshift 
dependent bias (which may be removed by using an rms 
effective redshift for each bin) and a scatter of ~3% to 
our L/ij estimates. These errors are small especially com- 
pared to systematics; therefore our samples of redshift 
binned galaxies are well represented by an average flux 
SED at the median redshift. 

4. RESULTS 

4.1. IR SFR Estimates 

SFR//? values that are obtained from CE01 template 
fits are illustrated in Figure 21 and tabulated in Table 
21 All of our redshift bins contain significant submillime- 
ter detections, which help to constrain the dust emis- 
sion peaks. We adopt CE01 template fits in the range 
A > 24/im for our preferred L/^ values. Though not for- 
mally the best x 2 , they are comparable to the best fits 
and including the 24 /xm photometry makes maximum 
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use of the available data. valu es obtained from Ta - 
ble [4] are consistent with results of Da ddi et al.l (|2005f >. 
who select BzKs to the same depth in the K band as 
presented here, and use MIPS 24 (im photometry to es- 
timate Ljr ~ 1.7 x 10 12 erg s" 1 for BzKs in the range 
1.4<z<2.5. 

4.2. UV SFR Estimates 

Table[5]illustrates the unweighted averages of estimates 
of UV SFRs from analysis of individual sBzK galaxies in 
each redshift bin. Averages excluded galaxies with poor 
fits to the spectral index, /3, identified by large x 2 (xi > 
2) or best fit /3 values that were pinned at the extreme 
of the allowed parameter range. In our three redshift 
bins 1.2 < z < 1.5, 1.5 < z < 2.0, and 2.0 < z < 3.3, 
these poor fit criteria excluded 41, 82, and 48 galaxies, 
respectively, from the averages. Average SFR^y r is in 
the range 12-285 M yr -1 , increasing with redshift. SFR 
in the highest redshift bin is affected by outliers; the 
median SFR^y r for individual galaxies are 10, 33, and 
106 M© yr -1 in each redshift bin, respectively. However, 
in keeping with the literature, we adopt the average of 
individual fits as our preferred indicator of UV SFRs for 
our sample. 

We also compute the SFRs from a single unweighted 
average spectrum of galaxies within each redshift bin. 
These estimates are systematically lower than the aver- 
ages of individual galaxies in each bin because the best 
fit UV continuum slopes to the average spectra indicate 
a lower dust correction than the average of individual 
fits. For redshift bins 1.2 < z < 1.5, 1.5 < z < 2.0, and 
2.0 < z < 3.3, the fits to average spectra had reduced 
X 2 values of 0.3, 1.8, and 3.8 respectively. We interpret 
these values to mean acceptable fits for the lower two 
redshift bins. 

There are eleven galaxies in the highest redshift bin 
with SFR >1000 M yr -1 . Checking the positions of 
these gala xies against the published LESS catalog f Weiss 
et al. 120091 ) indicates that thev are not submillimeter 
sources; separations between these galaxies and their 
nearest neighbor in the submillimeter catalog are all 
greater than 50". Five of these sources are detected 
in 24 /jm waveband, and their inferred luminosities and 
SFRs (from CE01 fits) are also high (SFR 24Aim >600 
M Q yr -1 ); however, as discussed below, SFR 2 4 Aim is 
known to be overestimated in this redshift and lumi- 
nosity range. One of these sources is detected in radio, 
with Li AG hz = 3 x 10 24 W Hz" 1 (SFRff^ z = 1650 
Mq yr -1 ), and therefore may be an example of previ- 
ously reported Optically Faint Radio G alaxies (OFRGs; 
IChapman et al.ll2004l : iCasev et aL"ll2009ft . 

IRAC colors are available for three of these eleven 
galaxies, and none of t hem a ppear in the AGN selection 
region of iStern et all (|2005T ) in IRAC color-color space. 
None of these sources are individually detected in X rays, 
although three of them are within the CDF-S (between 
5-10' from center). AGN contamination cannot be ruled 
out; however, we would expect obscuration of the AGN 
in rest-frame UV and optical wavebands to mean that 
star formation would dominate the emission (as opposed 
to the case in X-ray wavebands, where obscured AGN are 
a dominant confounding factor). These outliers may sug- 
gest either different dust physical properties or geometry 



Belli (|2003ft . Our estimates from Table ED 



in these galaxies. 

In c omparison with other literature works. [Daddi et aT] 
(2004) determine SFRs for a K Vega < 20 sample of 24 
sBzK galaxies in the GOODS-S field at z > 1.4 using 
SEP fitting and d ust correction using the method of 
IMeurer eT al. (1999), and find dust-corrected SFR in the 
range 100-600 M p yr -1 . In a spec t roscop ically selected 
sample of BzKs, lYoshikawa et all (|2010T ) find SFRs to 
vary widely, over three orders of magnitude. 

4.3. Radio SFR Estimates 

Radio fluxes, luminosities and associated SFRs are re- 
ported in Table ED and radio SFRs are compared to cal- 
ibrations in other wavebands in Table [8] It has been 
reported previously that SFR^g^ g exceeds SFRf^U, 
by approximately a factor of two (|BelIl I2003T) ■ D i screp - 
anci es betwe e n the radio SFR calibrations of iBell (2003) 
and ICondonl (| 19921 ) are not entirely surprising given the 
differ ent assumption s of each calibration. 

Daddi et al. (|2005f ) report radio stacking (weighted av- 
erage stacking + individual detections) of their Ky ega < 
20 sample to obtain a luminosity of 3.6 x 10 23 W Hz -1 , 
correspondin g to SFR ~ 210M Q yr -1 , using the radio 
calibration of Yun et aT] ([20011 ) . which is similar to our 
calibration of 

are consistent with these results. 

In stacking a Kab < 23 sample of BzKs, iDunne et al.l 
( 2009) reported a median sBzK luminosity of 1.28 xlO 23 
W Hz -1 corresponding to SFR=154±7 M© yr -1 , which 
is similar to our 1.2< z <1.5 bin result of 138±11 M Q 
yr , although a formal comparison is not possible be- 
cause of the 1.2 magnitude shallower depth of this present 
sample. Likewise, the values presented here are similar 
to results from the COSMOS survey (K s < 23 selected 
samp le; SFR in the range 30-100 M Q yr -1 IPannella et aT 
2009), where the radio-SFR calibration of lYun et al 
(|2001l ) is used. 

4.4. X-Ray SFR Estimates 

Our Chandra soft band stacked X-ray fluxes are in the 
range 3.2-9.5 xlO -18 erg s -1 cm -2 , see Table [7] Using 
observed frame, soft band fluxes and our assumed T — 
1.2 spectrum to convert flux to luminosity leads to rest 
frame, 2-10 keV luminosities in the range 7 — 50 x 10 40 
erg s -1 , indicated in Table [7] 

In comparison with other reported values of galaxies 
detected to t he sa me K band depth as presented here, 
IDaddi eta l. (2004) find rest frame, 2-10 keV luminosity 
of 8.6 x 10 41 erg s -1 (they use V = 2.1 in their flux 
to luminosit y conversion) in stac king 23 sBzKs in the 
K20 Survey dCimatti e t al. 2002J) that includes part of 
the CDF-S. IDaddi et all (|2005l ) find a rest frame, 2-10 
keV luminosity of 3.4 x 10 41 erg s -1 in stacking X-ray 
undetected sBzKs in 2 Ms Chandra data in GOODS-N 
(they use r = 2.0 in their flux to luminosity conversion). 
Thus we conclude that our stacked X-ray luminosities 
are consistent with other z ~ 2 star-forming galaxies 
reported in the literature. 

Our SFRs estimated from the rest frame 2- 1 keV lu- 
minosities, and the calibration of Ranalli ct al. (2003|) are 
in the range 15-100 M Q yr -1 , a fa ctor of » 5 low e r than 
the corresponding calib rations of iLehmer et al.l (2010) 
and iPersic et al.l (|2004D . As discussed in Section 13.41 
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SFR, 



Persic 
2-10fceV 



may be considered to provide upper limits. 



In the SFR calibration of lLehmer et al.l (|2010f ) , the stellar 
mass term contributes <23%, 7%, 3% and 1.4% to the 
SFRs in each redshift range of Table [SJ respectively. The 
trend of decreasing contribution from LMXBs to the to- 
tal X-ray luminosity as redshift increases is broadly con- 
sistent with models of the LMXB population and star 
formation histor y that predict the LM XB population to 
decline at z > 1 (IGhosh fe White! l200lh . 

5. DISCUSSION 
5.1. Comparison of SFR Estimates 

SFR estimates from X-ray through radio calibrations 
are compared in Table [5] The SFR/^ values in Table [3] 
are obtained from L] r estimates from fits of MIR-radio 
(A > 24/im) photometry to CE01 templates. To gauge 
their consistency, we plot the ratio of SFR computed 
from each calibration to SFR/# + ;yy in Figure [51 In this 
figure, SFR/i^/SFR/^+t/y is most nearly equal to one, 
reflecting that IR luminosity accounts for >90% of the 
total SFR in these galaxies. 

Comparing SFR 2 4 Aim to SFR/^+yy in Tableland Fig- 
ure [6] illustrates the overestimate of SFR at high redshift 
from this single waveband estimate. The resulting poor 
fit to the data in the highest redshift bin, 2.0 < z < 3.2, 
illustrated in Figure 21 overestimates L/^ by approxi- 
mately a factor of 6. Althoug h 24 /zm estimate s of Lj# 
can be robust at z < 1.5 (e.g. lElbaz et al.ll2010T ). overes- 
timation of Lir, particul arly at higher redshift has been 
previ ously reported (e.g. IPapovich et al.l 120071; Mur phy 
et al. 120091: iMuzzin et alJl2010t iNordon et al.N2010l: El- 
baz et al.l2011f). Using SvitzerIRS spectroscopy. Murphy 
et al. (|2009D concluded that this luminosity over-estimate 
arises due to unusually large polycyclic aromatic hydro- 
carbon features in these galaxies, and to a lesser extent, 
AGN contamination. 

Dust-corrected UV SFRs, shown in Table [5l from indi- 
vidual fits to sBzK photometry agree with SFRj#_|_[/y m 
the middle redshift bin, although the highest redshift bin 
exhibits larger SFR^y r and also larger error than lower 
redshift bins. This large SFR value is strongly affected 
by a small number of outliers with very high computed 
SFRs, which may suggest a modified extinction law for at 
least some galaxies above z ~ 2. Also, we find SFR^y r 
to underestimate SFR/ij + (yy in the lowest redshift bin 
by a factor of five. SFR^y r values were computed us- 
ing fits to the entire available broadband and medium 
band photometry (typically 12-14 sampled points in each 
fit); however, fits that were computed based upon broad- 
band photometry only (typically only 3 measured points 
in each fit) produced systematically higher SFRs. The 
medium band photometry clearly better samples the ob- 
served SED, and these results suggest a discrepancy be- 
tween SFR^y r and SFR/ fl+; yy for BzKs in this redshift 
range. 

Agreement between dust-corrected UV SFR and 
SFRjfl + j/y in BzKs to within a factor of ~ 2 has been 
report e d previously (e.g. iReddv et al.ll200 6; Dadd i et"aH 
120071* INordon et a l. 2010). IReddv et al.l 120101 ) find 



SFR^y r to agree with SFR determined fr om Ha spec 



troscopy for LBGs at z ~ 2. In particular. IReddv et al.l 
( 2006) fin d LBGs and BzKs w ith ages > 100 Myr to fol- 
low the iMeurer et al.l (|1999l ) relation while LBGs and 



BzKs with ages < 100 Myr have rest UV colors that are 
redder than expected for a given Lm + uy 

We have tried two methods of estimating radio fluxes: 
weighted average stacking of non-detections averaged 
with individual detections and median stacking of all 
sources. We find results of median stacking to be 
more consistent with SFR/^+^y, despite the fact that 
SFRjfl_|_j/y estimates were obtained with the method of 
weighted average stacking. Radio outliers in our sample 
may be due to residual AGN contamination. 

Among radio based SFR estimates, Figure [6] illustrates 
agreement between SFRf 1 1 q Hz and SFR/^+j/y to within 
a factor of two over our redshift range; this agreement 
is a consequence of the IR- Radio correlation for sBzKs. 
which is assumed i n the calibration of SFRf|g ff2 . The 
model of iCondonl (|1992l ) does estimate star formation 
from radio luminosity independent of the IR- Radio corre- 
lation; SFRf 4Q^™ exceeds SFR/^+^y by a factor of two 
over the observed redshift range, but the ratio of these 
two calibrations appears relatively insensitive to redshift. 
The discre panc ies between the radio SFR calibrations of 
iBel (|2003l ) and lCondonl (| 19921 ) are not entirely surprising 
given the different assumptions of each calibration. 

X-ray SFR indicators show wide variation in their 
ratios to SFR/ft+t/y. As discussed in Section 
14. 4i SFR^f^ l fe c e y may be interpreted as an upper 
limit. SFR^f i^ e e y yields estimates that are similar to 
SFR^iQ^y because of the subdominant contribution 
of the stellar mass term to the X-ray luminosity in 
these galaxies. Meanwhile SFRf°" a ^ y actually agrees 
with SFR/j^+yy to within a factor of two for z > 1.5. 
SFRg 'l^kev has been app lied to BzKs (e.g. iDaddi et all 
l2007al: lReddv et al.ll2005D. BX/BM galaxies (e.s. Reddv 
k Yy ml2QQ4|), and LBGs at z ~ 2 (e.g. |Reddy_et_aL| 
120101 ). In these studies, X-ray SFRs often agree with 
other waveband estimators, typically to within the same 
factors as reported here. 

However, interpretation of these X-ray SFRs depends 
upon an uncertain contamination fraction from obscured 
AGN. We speculate that AGN contamination may be 
present in our sample, particularly in the highest red- 
shift bin, in which we find L2-io/ceV ~ 10 erg s _1 . 
AGN contamination even at the level of 10% can re- 
quire d ownward adjustmen t to X-ray SFRs by a factor 
of 2-5 (jLehmer et al.ll2008D. Consequen t ly, X- ray SFRs 
would be over-estimated. iLehmer et aD (2008) compute 
a luminosity dependent AGN fraction in order to cor- 
rect X-ray stacking results in z ~ 3 galaxies; they find 
that ps 50 — 70% of the stacked 0.5-2 keV counts may 
arise from obscured AGN. If a similar AGN fraction ex- 
ists in our sample, then the X-ray SFRs would need to 
be adjusted downwards by a factor of ^2.5, bringing 
SFRi-iffkeV and SFR 2-iofceV into better agreement with 
other waveband indicators, and taking SFR^-ioleV out 
of agreement with other waveband indicators. 

5.2. Sources of Uncertainty 

In computing uncertainties in these SFR estimates, we 
have considered the effects of errors in photometry, spec- 
tral shape and redshift on luminosity estimates. These 
errors are reported in Table [8] and indicate a wide range 
in precision of the various waveband indicators. IR esti- 



Panchromatic Estimation of SFRs 



11 



mates are the most precise; X-ray SFR calibrations are 
the least precise because they rely upon empirical corre- 
lations with IR luminosity and thus introduce additional 
scatter into the SFR estimate. Each of these SFR cal- 
ibrations assume continuous star formation of at least 
10 s Myr and solar metallicity; therefore different assump- 
tions about timescales or chemical evolution cannot ac- 
count for systematic differences between the calibrations. 

We have excluded the uncertainty due to 
luminosity— SFR calibration. We are not able to assess 
the absolute uncertainties in these SFR calibrations be- 
cause in many cases they contain implicit dependencies 
on L IR -SFR calibration. SFR IR , SFR 24Mm , SFRf|^ ffz , 
SFRf^* v and SFR£f^ c e y depend i mplicit ly' upon 
the L/ fl -SFR calibration of iKennicuttl (|1998bl ). which 
has a reported systematic uncertainty of about a 
factor of 2-3. Similarly, SFR^™^ depen ds upon the 
L IR -SFR calibration of IB ell et all (|2005[ ). which also 
has a reported systematic uncertainty of a factor of 2. 
Similarly, SFR^y r depends upon a model dependent 
UV luminosity - SFR calibration, for whi ch various 
published values may differ by a factor of 2 (Kcnnicutt 
119981 . In the comparisons discussed here, we assume 
that the various SFR calibrations are consistent with 
each other and do not evolve with redshift. We compare 
them against each other, and disagreement can provide 
evidence of systematic offsets in a given calibration. 

For generalizing the results of this sample of BzKs 
to the general population of BzKs, errors may be esti- 
mated from bootstrap resampling, which would incorpo- 
rate sample variance, and undoubtedly increase the size 
of random errors. We do not consider this sample vari- 
ance here, because our primary aim is to compare the 
SFR calibrations to each other, to assess their consis- 
tency, rather than comparing SFRs of BzKs as a popu- 
lation to other populations of star-forming galaxies. 

We consider the differing sensitivities to star formation 
among the various wavebands. We estimate the lowest 
average SFR detectable in each waveband and redshift 
bin from the reported image depths by computing the av- 
erage flux that would yield a 3a stacked detection given 
our sample sizes. Converting this flux to a minimum 
detectable average luminosity at each bin median red- 
shift yields SFR values < 3,~ 3, and ~ 15 MQyr -1 for 
SFRg" corr , SFR 24Mm , and SFR flodio respectively, at all 
sampled redshifts. SFR^y r and SFR/^ are even more 
sensitive than SFRyy corr and SFR24 Aim ; quantitative es- 
timate would require detailed modeling of the SED fitting 
procedure. The X-ray calibrations are less sensitive, with 
minimum detectable SFR in the range 30-66 and 150-330 
M Q yr _1 (increasing with redshift) for SFR^"g" e l y and 

SFR^f iQfcev respectively. X-ray SFR estimators have low 
sensitivity and are subject to bias due to residual AGN 
contamination and it is not clear which of these effects 
explains the overestimation of SFR. 

6. CONCLUSION 

The main results of this paper are summarized in the 
comparison of SFR indicators given by Table [5] and Fig- 
ure |51 We consider SFR determined from panchromatic 
estimation of L/^j the most comprehensive approach, 
where the available data exist, notwithstanding the chal- 
lenges of accurate fitting of spectral templates. With this 



method, average SFR of redshift binned galaxies can be 
determined with ^5%— 10% random uncertainty (though 
considerably larger for individual sources); however, the 
total error is dominated by a factor of ^2 systematic 
uncertainty. 

This systematic uncertainty, not included in the er- 
rors indicated here, of a factor of ~ 2 is present in 
the Li R — SFR and Ljjy— SFR calibrations that underlie 
each method of SFR estimation. Each of the calibra- 
tions discussed in this paper either implicitly or explic- 
itly assumes a continuous star formation history and so- 
lar metallicity; consequently differing assumptions about 
these parameters should not account for systematic dif- 
ferences. 

We find that dust-co rrected UV SFR, using the method 
of lMeurer et al.l (|1999f ). agrees with SFR/^+c/y for galax- 
ies in the range 1.2 < z < 1.5, but over-estimates 
SFR/fl: + [/y in the range 2.0 < z < 3.2 due to a small 
population o f outl i ers. Radio SFRs estimated from the 
calibration of lBelll (|2003f ) are in agreement with the total 
SFRs to within errors. The SFR calibration of iCondonl 
( 1992) overestimates the total SFR by a factor of two for 
BzK galaxies over this redshift range. 

Perh aps surprisingly, the calibration of IRanalli et al.l 
(2003) yields SFR estimates that agree with SFR/ij+yy 
and other indicators in the range 1. 5 < z < 3.2 , while 
the better suited X-ray c alibrations of lPersic et al.l (|2004f ) 
and lLehmer et all ([2010D over-estimate SFR in the range 
1.5 < z < 3.2. One possible explanation is mild 
AGN contamination in the upper two bins, which would 
bring SFR^f ™ c ey and SFR^o& into rou g h agreement 
with SFR/ fl+[ /y and cause SFR^u^ kiv to underestimate 
SFR. X-ray SFR estimates are also notably less sensitive 
than other waveband indicators. 

There is an evident variation in the accuracy and preci- 
sion of SFR calibrations, and each method has its limita- 
tions and caveats. Radio and X-ray SFR calibrations rely 
upon empirical correlations of flux in their wavebands to 
Lir that introduce inevitable scatter, particularly in the 
X-ray calibrations. 24 /im only and UV SFRs work for 
most galaxies but have important exceptions and cannot 
be applied universally. 

Our analysis of Kab < 21.8 sBzKs in the redshift range 
1.2 < z < 3.2 confirms that they are IR luminous, star- 
forming galaxies for which approximately 90% of the to- 
tal star formation is obscured by dust. By fitting to CE01 
templates, we find average IR luminosities for redshift 
binned sBzK galaxies at median redshifts 1.4, 1.8, 2.2 to 
be 3.0±0.3x 10 11 L , 4.0±0.4x 10 11 L Q , 8.3±0.7x 10 11 
L Q . We find SFRjij+uv at these redshifts to be 55±6, 
74±8, 154±17 M yr" 1 respectively. 
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7. APPENDIX: ERROR ESTIMATION IN STACKING 

7.1. Ordinary Average 

Section [3] presents the method used in this paper for 
estimating the aggregate flux, p, from N prior positions 
that include a combination of / individual detections, x^, 
and a stacking detection, xs, from S undetected sources, 
where N — I + S, given by 



(18) 



The variance of. of this estimate is related to the indi- 
vidual errors <7j and o~ s as well as their covariances, o~i ]S 
from standard error analysis 




i=l 
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dp 



\p„ 2 dfi dp 
^ " -"<);■, dx a 

I— 1 



(19) 

We assume that the covariances between individual de- 
tections and the stacking detection, o~i iS and the covari- 
ances between separate individual detections, cn,j, are 
zero. 
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(22) 



Equations[TH]and[22]are used to compute the average and 
error, respectively, of individual and stacking detections 
reported in this paper. 

In the case of all measurement errors being equal, 
which is a good approximation in the case of LESS data, 
then Oi = a and a s = a/y/N. Using Equation [22j is 
computed as 



1 

TV' 



S 2 - 



s 



I + S 



y/la 2 + Sa 2 



7.2. Weighted average: 



(23) 

(24) 
(25) 



An alternative approach to combining individual and 
stacking detections into a single aggregate flux estimate 
is to use a weighted average of individual and stacking 
detections. Assume there are N = I + 1 flux measure- 
ments consisting of / individual detections and a single 
stacking detection. Each individual flux measurement 
and the stacked flux measurement is considered as an 
independent flux measurement for the purpose of com- 
puting an average, and these measurements are combined 
as an inverse variance weighted average. 



N 



v-JV i 



(26) 



(27) 



Equations 1261 and 1271 are also used to estimate the aggre- 
gate flux and error for representative data reported in 
this paper. In particular, the errors computed according 
to this method are numerically equal to the errors com- 
puted from the ordinary average, Equation 1221 to better 
than three significant digits. 

For clarity, Equation [57] can be written to include 
stacking and individual detections separately 



1 



(28) 



In the case of all measurement errors being equal, Equa- 
tion [28] can be simplified using <jj = a and a s — a/y/N, 
just as in the ordinary average error computation: 



v / 1 



i + S 



y/T+S 



(30) 
(31) 



Thus for the case of all measurement errors being equal, 
the error of the weighted average is identical to the error 
of the ordinary average. 

7.3. Comparison of ordinary and weighted average: 

In the case of identical errors for the individual mea- 
surements, o~n = cr, then weighted and ordinary averages 
give the same result for the aggregate flux. In the case 
of measurement errors being unequal, then the weighted 
average will in principle have the smaller error; however, 
the difference will be small, and if the errors are not in- 
dependent, e.g. if brighter sources have larger errors, 
then the weighted average introduces a bias to the flux 
estimate p' . For instance, if dim sources are always mea- 
sured with better precision than bright sources, then a 
weighted average of the population of all sources will al- 
ways be biased toward dim sources. This circumstance 
could arise if flux measurement errors are dominated by 
Poisson counting statistics. However, if the flux measure- 
ment errors are uncorrelated with the flux, then there will 
be no problem with the weighted average. 
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Fig. 1. — BzK diagram illustrating corrected NIR-optical colors (z — K)ab vs. (B — z)ab f° r MUSYC K < 21.8 selected sources (gray 
points), including redshift binned, star-forming galaxies (with X-ray detected AGN removed, see the text) in the range 0.9 < z < 1.2 (red 
circles), 1.2 < z < 1.5 (green triangles), 1.5 < z < 2.0 (blue squares) and 2.0 < z < 3.2 (purple stars). The sBzK region is located above 
the diagonal line, and is defined as BzK > 0.2. The pBzK region is the wedge shaped area above the horizontal dotted line, z — K > 2.5. 
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Fig. 2. — Histograms of photometric and spectroscopic redshifts for sBzK galaxies in ECDF-S. The upper histogram indicates photometric 
redshifts. The lower filled (red), histogram indicates spectroscopic redshifts. The vertical dotted line demarcates z < 1.2 galaxies that are 
excluded from analysis. 
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Fig. 4. — Spectral energy distributions of redshift binned sBzK galaxies obtained from stacking analysis, and fits to various models. 
Spectral flux density, in units of /iJy, is plotted vs. observed frame wavelength, in units of Angstroms. UV through radio flux measurements 
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dashed (blue) CE01 template fitted to the A > 24/im spectrum, and solid (red) CE01 template fitted to only the 24 (im data point. For 
each model fit, the region of the spectrum used for the fit is indicated by the solid portion of the curve. Top: galaxies in the redshift range 
1.2 < 2 < 1.5; middle 1.5 < z < 2.0; bottom 2.0 < z < 3.2. See also TableHfor fit details. 
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Fig. 5. — Left panel illustrates simulated, observed frame spectra for sBzK galaxies in the range 1.5 < z < 2.0. The dotted curve indicates 
the best fit CE01 template placed at the median redshift, z = 1.756. The solid curve indicates the result of averaging 233 identical CE01 
template sources, distributed in redshift according to the observed distribution of sBzKs, with photometric redshift errors added according 
to the distribution shown in Figure [3] The right panel illustrates the frequency distribution of hm fractional errors due to photometric 
redshift errors and dispersion from individual galaxy SEDs, for sBzK galaxies in the range 1.5 < z < 2.0. Solid curve (red) indicates a 
Gaussian fit with mean = 0.14 and tr=0.06. Thus Lj^ for this redshift bin is overestimated by 14%; photometric redshift and galaxy SED 
dispersion contribute a scatter of 6%. 
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TABLE 1 

Summary of Results of AGN Selection 
Criteria Applied to sBzK Sample. 



Test 


Detected 


AGN 


X-ray 




IRAC 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


X-ray 


110 


107 


83 


7 


18 


1IR 


35 


22 


7 


15 


2 


IRAC 


649 


25 


18 


2 


7 



Note. — Column (1) indicates the waveband 
test used for AGN discrimination. Column (2) in- 
dicates the number of sources with data in each 
respective waveband. Column (3) indicates the 
number of AGN confirmed by each waveband test 
irrespective of tests in other wavebands. Entries 
in Columns (4-6) indicate the number of sources 
classified as AGN according to each waveband as 
follows: Diagonal entries (in boldface) indicate 
sources that are uniquely identified as AGN in 
only one waveband. Off diagonal entries indicate 
sources that are identified as AGN in at least the 
two wavebands indicated by the respective row and 
column headings. 
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TABLE 2 

Redshift Binned Sample of Non-AGN Star-Forming Galaxies (sBzKs) in the ECDF-S 



Total Sample Individual Detections 



Redshift 
Range 
(1) 


Median 
Redshift 
(2) 


Number 

( z phot> Z s p ec ) 

(3) 


MIPS 
24 /im 
(4) 


MIPS 
70 £tm 

(5) 


LESS 
870 Atm 

(6) 


VLA 
1.4 GHz 

(7) 


Chandra 
0.5-2 keV 
(8) 


1.2 < z < 1.5 


1.383 


156 (152,4) 


64 


4 


2 


2 


1 


1.5 < z < 2.0 


1.753 


215 (202,13) 


90 


9 


6 


9 


1 


2.0 < z < 3.2 


2.272 


139 (122,17) 


61 


1 


1 


5 






Note. — Redshift binning scheme selected for this analysis. Column (1) indicates the redshift 
range for each bin. Column (2) indicates the median redshift of the sampled galaxies in each 
bin. Column (3) shows the total number of galaxies in each bin, with numbers of photometric 
and spectroscopic rcdshifts, respectively, in parentheses. Columns (4-8) indicate the number of 
sources that are individually detected in each waveband; fluxes from individual detections are 
combined with stacked fluxes of the remaining sources in each redshift bin, as discussed in the 
text and Appendix. 
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TABLE 3 

UV-Radio Average Flux Densities for Redshift Binned sBzKs 







1.2 < 2 < 1.5 


1.5 < z < 2.0 


2.0 < 


z < 3.2 


Band 


A 


S„ 


a S 




a S 


S„ 


°"S 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


T T 

u 


U.oO 


0.45 


0.06 


0.36 


0.05 




U.UO 


JJ 


0.46 


0.61 


0.04 


0.58 


0.04 


0.59 


0.04 


V 


0.54 


0.65 


0.04 


0.67 


0.04 


0.73 


0.05 


H 


0.65 


0.85 


0.05 


0.80 


0.05 


0.86 


0.05 


I 


0.86 


1.57 


0.17 


1.36 


0.16 


1.16 


0.16 


z 


0.90 


1.63 


0.24 


1.34 


0.23 


1.18 


0.24 


J 


1.25 


3.77 


0.63 


4.31 


0.63 


2.79 


0.61 


H 


1.65 


4.25 


0.79 


4.97 


0.80 


3.52 


0.71 


K 


2.13 


9.70 


1.24 


10.32 


1.22 


9.69 


1.24 


24 fim 


24 


94 


2 


110 


2.0 


120 


2.1 


70 /im 


70 


320 


81 


340 


83 


440 


85 


250 /im 


250 


2000 


880 


9080 


750 


1400 


930 


350 (im 


350 


2020 


690 


5700 


590 


920 


730 


500 Mm 


500 


1030 


480 


4300 


400 


1300 


510 


870 fim 


870 


240 


94 


510 


80 


530 


100 


1.4 GHz 


214000 


15 


0.8 


12 


0.7 


15 


0.8 


610 MHz 


490000 


27 


5.2 


32 


4.3 


15 


5.3 



Note. — Column (1) — waveband. Column (2) — effective, observed 
frame wavelength in units of (im. Columns (3,5,7,9) — average, observed 
flux density in units of £tjy for each redshift bin. Column (4,6,8,10) — error 
in flux density in units of /ijy for each redshift bin. 
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TABLE 4 

Redshift Binned sBzK IR-Radio Fit Summary 



Redshift 
(1) 


Fit Type 
(2) 


(3) 


X 2 (df) 
(4) 


SFR/fl 
(5) 


1.2 < z < 1.5 


CEOl (>24 fim) 
CEOl (>24 /tin) 
CEOl (24 /im) 


3.0±0.3 
1.9±0.2 
4.1±0.3 


48.23(7) 

5.92(6) 

3.28(0) 


51±5 
33±3 
70±9 


1.5 < z < 2.0 


CEOl (>24 fim) 
CEOl (>24 /mi) 
CEOl (24 /im) 


4.0±0.4 
5.0±0.9 
8.6±0.5 


100.94(7) 

75.31(6) 
0.03(0) 


68±7 

86±15 

148±45 


2.0 < z < 3.2 


CEOl (>24 /mi) 
CEOl (>24 /mi) 
CEOl (24 (im) 


8.3±o.r 

6.3±0.7 
35.4±2.1 


51.24(7) 
13.23(6) 
1.82(0) 


143±12 
108±12 
608±91 



Note. — Column (1) indicates redshift range for the sample. 
Column (2) specifies the type of model fit (see text). Column (3) 
indicates the computed Lj/j from the best fit model, in units of 10 
Lq . Column (4) indicates the \ value for the best fit, with degrees 
of freedom in parentheses. Column (5) indicates the SFR, in units 
of Mq yr -1 , computed from the fit-derived IR luminosity. Random 
errors correspond to 68% confidence intervals, and do not include 
substantial systematic error, as discussed in the text. 
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TABLE 5 

Redshift Binned sBzK UV Luminosity and SFR Estimates 







Average of Individual Spectra 




Average Spectrum 


Rcdshift 
(1) 


Pfit 

(2) 


log(IRX) 
(3) 


Ai600 
(4) 


SFR^" corr 
(5) 


SFR<^ rr 
(6) 


SFRg^T r 
(7) 


1.2 < z < 1.5 
1.5 < 2 < 2.0 
2.0 < 2 < 3.2 


-1.50 
-0.80 
-0.59 


1.12 
2.15 
2.21 


1.44 
2.84 
3.26 


4±0.1 
6±0.3 
11±0.8 


12±1 
64±4 
285±30 


10±3 
36±12 
65±28 



Note. — Column (1) indicates sample redshift range. Columns (2-6) correspond to 
unweighted averages of results from fitting individual galaxies in each redshift bin. Column 
(2) indicates the best fit slope, /?, to the fx spectrum where fxocX 13 . Column (3) indicates 
the IR-UV ratio, log(S/ij/Si6oo)- Column (4) indicates the attenuation in magnitudes at 



1600 A derived from the best fit spectral index, 0. Column (5) and (6) refer to SFRs in 
units of Mq yr -1 . Column (7) indicates the SFR from IRX-/3 correction applied to the 
corresponding (unweighted) redshift bin averaged spectrum. 
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TABLE 6 

Redshift Binned sBzK Radio Flux, Luminosity and SFR Estimates 



Redshift 


Si. 4 


So. 610 


Li. 4 x 10 22 




SFRf|" 


Range 


A»Jy 




W Hz" 1 


M yr" 1 


M yr" 1 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1.2<z<1.5 


14.6 (12.1)±0.8 


26.6 (14.8)±5.2 


13.9 (11.5)±0.8 


166 (138)±11 


76 (64)±5 


1.5<z<2.0 


11.7 (8.9)± 0.7 


32.2 (33.2)±4.3 


19.4 (14.7)±1.1 


232 (176)±16 


107 (81)±7 


2.0<z<3.3 


15.4 (8.4)± 0.8 


14.6 (17.2)±5.3 


46.6 (25.3)±2.5 


560 (303)±36 


257 (139)±17 



Note. — Column (1) indicates the sample bin redshift range. Columns (2) and (3) include 
individual detections and weighted average stacked flux at 1.4 GHz and 610 MHz respectively. In 
Columns (2,3,4-6), results of median stack of all sources are indicated in parentheses. Column 
(4) indicates the derived re st frame 1.4 GHz luminosity. Column (5) indicates the computed SFR 
according to Condon (1992). Column (6) indicates the computed SFR according to Bell (2003). 
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TABLE 7 

Redshift Binned sBzK X-ray Flux, Luminosity and SFR Estimates 



Redshift 
Range 
(1) 


Flux 
0.5-2 keV 
(2) 


Flux 
2-8 kcV 

(3) 


Luminosity 
Rest 2-10 keV 
(4) 


C"P"D Ranalli 
orn 2 _ 10feeV 

M Q yr" 1 
(5) 


CTjn Persic 
arrv 2-10fceV 
Mq yr" 1 
(6) 


QT?yj Lehmer 
or rv, 2 -10keV 
Mq yr" 1 
(7) 


1.2 < z < 1.5 


3.2±1.0 


31±21 


7±7 


15±22 


73±67 


62±46 


1.5 < 2 < 2.0 


7.4±1.2 


44±16 


29±10 


58±75 


292±116 


291±123 


2.0 < z < 3.2 


9.5±1.2 


27± 6 


50±14 


100±126 


499±170 


507±203 



Note. — Column (1) — the sample bin redshift range. Column (2) — observed flux density in 
Chandra soft band (0.5-2 keV), in units of 10~ 18 erg s~ 1 cm -2 . Column (3) — observed flux density 
in Chandra hard band (2-8 keV), in units of 10 -18 erg s - 1 cm -2 Column (4) — the luminosity in rest 
frame (2-10 keV), computed from observed, soft band flux, assuming spectrum with photon index, 
T=1.2 and E e =20 keV, at the med ian redshift. Units are 10 40 erg s _1 . Colu mn (5) - SFR from 
the method of IRanalli et all f2003|) . Column (6) - SFR from the method of IPersic et~aT] p004fi . 
Column (7) - SFR from the method of lLehmer eTaLl (2010). All SFRs in units of M yr 1 assuming 
Salpeter IMF (0.1 - 100 M ). 



TABLE 8 

Redshift Binned sBzK Star Formation Rate Summary 



Redshift 
(1) 


SFR IR 
(2) 


SFR24pm 

(3) 


SFRg^ r 
(4) 


Condon 
artt 1.4GHz 

(5) 


(6) 


c T7 1 'p Ranalli 
or ^ 2 -10fceV r 

(7) 


CT^ry Per sic 
or ^2-10fceV 

(8) 


cppLe/imer 
or ^2-10feeV 

(9) 


SFR/ij + (7v 
(10) 


1.2 < z < 1.5 


51±5 


79±9 


12±1 


138±11 


64±5 


15±22 


73±67 


62±46 


55±6 


1.5 < z < 2.0 


68±7 


172±51 


64±4 


176±16 


81±7 


58±75 


292±116 


291±123 


74±8 


2.0 < z < 3.2 


143±12 


620±85 


285±30 


303±36 


139±17 


100±126 


499±170 


507±203 


154±17 



Note. — Star formation rates for redshift binned sBzKs. Column (1) — sample redshift range. Column (2) — SFR, from integrated IR 
luminosity using CE01 template fits to MIR-radio photometry (A > 24/jm). Column (3) — SFR from CE01 template fit to 24 (im data 
only. Column (4) — SFR from average of i ndividua l fits t o UV cont i nuum with dust correction. Columns (5) and (6) — SFR from median, 
stacked 1.4 GHz flux using the method o flUo ndonl 0-9921) a nd [Bell (2003) respectivel y. Columns (7-9) — SFR from soft band X-ray data 
using the method of Ranalli ct al. (2003), Persic ct al. (2004) and Lchmcr ct al. (2010) respectively. Column (10) — the total SFR from the 
sum of Column (2) and corresponding uncorrected UV SFRs taken from Table [5] Units are Mq yr — 1 for all columns. Values in boldface 
fall within 2<r of the best estimate, SFR/H+yv"- Errors do not include substantial systematic uncertainty, as discussed in the text. 



